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Low-lying states in nuclei are investigated using an ensem- 
ble of random interactions. Both in the nuclear shell model 
and in the interacting boson model we find a dominance of 
J p = + ground states. It is shown that this feature is not 
due to time reversal symmetry. In the shell model, evidence 
is found for the occurrence of pairing properties, and in the 
interacting boson model for both vibrational and rotational 
band structures. Our results suggest that these features rep- 
resent general and robust properties of the model space, and 
do not depend on details of the interactions. 

PACS numbers: 21.60.Cs, 21.60.Ev, 21.60.Fw, 24.60.Lz 



I. INTRODUCTION 

Despite the fact that nuclei are complex many-body 
systems with many degrees of freedom, their spectral 
properties often show very regular features. A recent 
analysis of experimental energy systematics of medium 
and heavy even-even nuclei suggests a tripartite classifi- 
cation of nuclear structure into seniority, anharmonic vi- 
brator and rotor regions Jl],|| . Plots of the excitation en- 
ergies of the yrast states with J p = 4 + against J p = 2 + 
show a characteristic slope for each region: 1.00, 2.00 and 
3.33, respectively. In each of these three regimes, the 
energy systematics is extremely robust. Moreover, the 
transitions between different regions occur very rapidly, 
typically with the addition or removal of only one or two 
pairs of nucleons. The transition between the seniority 
region (either semimagic or nearly semimagic nuclei) and 
the anharmonic vibrator regime (either vibrational or 7 
soft nuclei) was addressed in a simple schematic shell 
model calculation and attributed to the proton-neutron 
interaction || . The empirical characteristics of the col- 
lective regime which consists of the anharmonic vibrator 
and the rotor regions, as well as the transition between 
them, have been studied ||,|5| in the framework of the 
interacting boson model (IBM) || . An analysis of phase 
transitions in the IBM [|],[| has shown that the collec- 
tive region is characterized by two basic phases (spheri- 
cal and deformed) with a sharp transition region P,^fj[, 
rather than a gradual softening which is traditionally as- 
sociated with the onset of deformation in nuclei. 

In a separate development, the characteristics of low- 
energy spectra of many-body even-even nuclear systems 



have been studied recently in the context of the nuclear 
shell model (SM) with random two-body interactions 
JTT] , p~2| | . Despite the random nature of the interactions, 
the low-lying spectra still show surprisingly regular fea- 
tures, such as a predominance of J p — + ground states 
separated by an energy gap from the excited states. This 
is contrary to the traditional wisdom in which the favor- 
ing of J p = + ground states is attributed to the nu- 
clear pairing arising from the short-range nuclear force. 
A subsequent analysis of the pair transfer amplitudes has 
shown that pairing is a robust feature of the general two- 
body nature of shell model interactions and the structure 
of the model space Jjj|. On the other hand, no evidence 
was found for rotational band structures. 

The existence of robust features in the low-lying spec- 
tra of medium and heavy even-even nuclei [[^] suggests 
an underlying simplicity of low-energy nuclear structure 
never before appreciated. In order to address this point 
we carry out a study of the systematics of energy levels 
in the framework of the SM and the IBM with random 
interactions. We address time-reversal symmetry in con- 
nection with the dominance of J p — + ground states, 
look for regular spectral properties, and investigate the 
effect of many-body interactions. 



II. THE NUCLEAR SHELL MODEL 

We first consider the properties of nuclei in the pres- 
ence of random interactions within the context of the 
shell model. As a model space we take that of N identical 
nucleons in the sd shell, which consists of single-particle 
orbitals with j = 1/2, 3/2 and 5/2 . The case of A r = 6 
particles is one of the examples considered in [ pdp2] ] and 
referred to as corresponding to the nucleus 22 O. For iden- 
tical particles the isospin is the same for all states, and 
hence does not play a role. In Ref. E^] it was shown 
that the single-particle energies have little effect on the 
results, and therefore they are not considered here. The 
two-body interactions can be expressed as 
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The coefficients Cl ; correspond to the 30 two-body ma- 
trix elements for identical particles in the sd shell. They 
are chosen independently from a Gaussian distribution of 
random numbers with zero mean and variance v 2 . 



(3) 



Here <> denotes an ensemble average. The ensemble of 
Eq. (||) satisfies the requirement that it is invariant under 
orthogonal transformations (i.e. a change of basis). The 
variance of the Gaussian distribution v 2 is independent 
of the angular momentum and only represents an over- 
all energy scale. All results that we present here were 
obtained using 1000 runs. 

For N = 2 particles the Hamiltonian matrix is entirely 
random and the ensemble coincides with the Gaussian 
Orthogonal Ensemble (GOE), which is characterized by 
a semi-circular distribution of eigenvalues |l4| 
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whose width depends on the dimension d of the Hamil- 
tonian matrix |15| according to 
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In Table | we show the percentage of the total number of 
runs for which the ground state has a given angular mo- 
mentum. Clearly, the angular momentum for which the 
width of the distribution is largest will be the most likely 
to be the ground state. In this noted earlier, the 

widths depend directly on the corresponding dimension 
of the basis. Thus, the J — 2 state is the most likely 



to be the ground state, followed by J = and then by 
J = 1,3,4, exactly as seen in the table. These various 
points are made clearer from the semi-circular level dis- 
tributions shown in Fig. ^. 

For N > 2 particles the ensemble is the Two-Body 
Random Ensemble (TBRE) in which the ^-body 

matrix elements are correlated and can be expressed in 
terms of the random two-body matrix elements of Eq. (|^) 
by the usual reduction formulae. The eigenvalues now 
follow a Gaussian distribution [fill vM 
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As an example, we consider the nuclei 20 > 22 O with four 
and six valence neutrons. In Table [n] we show the per- 
centage of the total number of runs for which the ground 
state has a given angular momentum. The corresponding 
level distributions are shown in Figs. || and ||. For N = 4 
particles the percentage of J p = + ground states is 55.9 
%, significantly larger than for N — 2 (see Table Q). Note, 
however, that the percentage of + states in the model 
space is only 11.1 % for this case. Similar results hold for 
TV — 6 particles (see the last column of Table ||) ; namely 
the percentage of + ground states dramatically exceeds 
the percentage of + states in the basis. Our results for 
N = 6 arc in agreement with those obtained earlier by 
Johnson et al. |yjl§. 

In the cases of N = 4 and 6, the percentage of ground 
states associated with each angular momentum is also 
correlated with the widths of the distributions, which 
are now Gaussian (see Figs. || and ||). The key difference 
between these results and those for N = 2 is that here 
there is no direct connection between the width and the 
size of the basis for a given angular momentum. J p = + 
ground states predominate for N = 4 and 6 even though 
they have much fewer basis states than some of the other 
angular momenta. 

The observed preponderance of J p = + ground states 
for N > 2 is surprising, considering that there is no ob- 
vious pairing character in the assumed random forces. 
Thus the question remains: what is it that produces 
this dominance of + ground states in even-even many- 
body systems? One possibility is that it arises because of 
the time-reversal invariance of the random Hamiltonian. 
Since time-reversed states play an important role in the 
formation of correlated + (Cooper) pairs which in turn 
can give rise to favored collective many-body states, it is 
conceivable that time-reversal invariance may contain a 
built-in preference for J p = + many-body ground states 

s 

To see whether this is indeed the case, we consider 
what happens when we break time-reversal invariance in 
the random two-body interactions. This can be done by 
introducing a Gaussian Unitary Ensemble (GUE) rather 
than a Gaussian Orthogonal Ensemble (GOE) to ran- 
domly generate the two-body matrix elements. More 
specifically, we consider a two-body Hamiltonian of the 
form Iilj20tl 
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The coefficients and tjl^ are chosen independently 
from a Gaussian distribution of random numbers with 
zero mean and variance v 2 as 

(Ci«Ct; v ) ^S w (i + S t]tl , r )v 2 , 

(pLij-HL'.,., ) = &LL' (1 - 8ii,i'j>)v 2 , 

(C^ v ) - 0. (8) 

For e = and 1 they correspond to GOE and GUE, re- 
spectively. The Hamiltonian is time-reversal invariant if 
the two-body matrix elements are real, i.e. e = 0. The 
breaking of time-reversal symmetry can be studied by 
taking < e < 1. For a given value of J, the above en- 
semble for two-body interactions gives a semicircle level 
density. The normalization is chosen such that the radius 
of this semicircle distribution does not depend on e [SI . 
In Table III we present the results for N = 6 identical 



particles in the sd shell. For < e < 1 the time-reversal 
invariance is broken. We see that the dominance of + 
ground states increases with e from 67.7 to 76.8 %. On 
the basis of these results, we conclude that time-reversal 
invariance of the two-body interactions is not the origin 
of the dominance of + ground states. 

For the cases with a J p = + ground state we calculate 
the probability distribution of the energy ratio 
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E(A+)-E(0+) 
E(2+)-E(0+) 
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As noted earlier, this energy ratio has characteristic val- 
ues of R « 1, 2 and 10/3 for the seniority, vibrational and 
rotational regions, respectively. In Fig. |we show the re- 
sults for N = 6 particles. The probability distribution 
shows a broad peak between 1 < R < 2, with a maxi- 
mum around 1.3. This suggests that on average a system 
of identical nucleons tends to behave in accord with the 
seniority regime of QJ^]. An analysis of the amplitudes 
for pair transfer between ground states has shown that 
pairing is a robust feature of two-body shell model inter- 
actions and arises from a much broader class of Hamilto- 
nians than the ones usually considered Jl3| | . This finding 
and ours are in qualitative agreement with the empiri- 
cal observation of very robust spectroscopic features in 
the seniority regime (l],^] . Since the distribution extends 
to R = 2, we conclude that there is also some evidence, 
although minimal, for vibrational structure in our cal- 
culations. On the other hand, there is no evidence for 
the occurrence of rotational bands, at least within the 
context of the model space we have considered (see also 
0)- 



In the next section we carry out a study of the system- 
atics of collective levels in the IBM with random interac- 
tions, and look for evidence for vibrational and rotational 
bands. 



III. THE INTERACTING BOSON MODEL 

In the IBM collective nuclei are described as a sys- 
tem of N interacting monopole and quadrupole bosons 
H . The one-body Hamiltonian of the model contains the 
single-boson energies 



Hi = £q ■ S + 62 eft ■ d 
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and the two-body Hamiltonian contains the various two- 
boson interactions 
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with 



= (,t x S t)(0)/V2 , 
P f 2 = (dt x rft)(o)/^2 , 

Pl = ( S t X dt)(2) , 

Pl = (^ x rft)(2) /V 2 , 
Pl = (d+ x dt)( 4 VV2 . 



(12) 



The coefficients e j, and fiy correspond to the 2 one-body 
and 7 two-body matrix elements, respectively. They are 
chosen independently from a Gaussian distribution of 
random numbers with zero mean and variance v 2 accord- 
ing to 



(££££') = S LL > 2v 2 , 

(di«Ctj ) = &lv (i + kit'?) v 2 
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First we consider the most general one- and two-body 
IBM Hamiltonian 
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Note that to remove the N dependence of the matrix 
elements of fc-body interactions, we have scaled Hj~ by 

rL=i(-^ + 1 — l n a tl calculations we take = 16 and 
1000 runs. For each set of randomly generated many- 
body matrix elements we calculate the entire energy spec- 
trum and the B(E2) values between the yrast states. 

Just as in the case of the nuclear shell model [[ll] , we 
find a predominance of J p — + ground states; 63.4 % 
of the ground states have this value of the angular mo- 
mentum, even though only 3.3 % of the basis states do. 
Other angular momenta for which there are relatively 
high ground-state probabilities are J p = 2+ (13.8 %) 
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and J p = 32 + - the maximum value of the angular mo- 
mentum (16.7 %). 

For those cases having a J p = + ground state we have 
calculated the probability distribution of the energy ratio 
R of Eq. (^). Fig. || shows a remarkable result: the prob- 
ability distribution P(R) has two very pronounced peaks, 
one at R ~ 1.95 and a narrower one at R ~ 3.35 
These values correspond almost exactly to the harmonic 
vibrator and rotor values of 2 and 10/3 (see Table IV). 

Energies by themselves are not sufficient to decide 
whether or not there exists a collective structure. Levels 
belonging to a collective band are connected by strong 
electromagnetic transitions. In Fig. |^ we show a cor- 
relation plot between the ratio of B(E2) values for the 
4^ — > 2+ and 2f — ► 0^ transitions and the energy ratio 
R. For the B{E2) values we use the quadrupole operator 



Q Al ( X ) = ( S td + dt s )(2) +x(d tJ ) (2) 



(15) 



with x — — v7 /2. For completeness, in Table IV we show 
the results for the three symmetry limits of the IBM |j) . 
In the large N limit, the ratio of B(E2) values is 2 for the 
harmonic vibrator and 10/7 both for the 7- unstable rotor 
and the pure rotor. There is a strong correlation between 
the first peak in the energy ratio and the vibrator value 
for the ratio of B{E2) values (the concentration of points 
in this region corresponds to about 50 % of all cases), and 
for the second peak and the rotor value (about 25 % of 



all cases) 21 



The results presented in Figs. || and || were obtained 
with random interactions, with no restriction on the sign 
nor the magnitude of the one- and two-body matrix el- 
ements. It is of interest to make a comparison with a 
calculation in which the parameters are restricted to the 
'physically' allowed region. To this end, we consider the 
consistent-Q formulation |E| which uses the same form 
for the quadrupole operator, Eq. (pf), i.e. with the same 
value of x f° r the E2 operator and the Hamiltonian 



H = eh d - kQ(x) ■ Q(x) ■ 



(16) 



The parameters e and n are restricted to be positive, 
whereas x can be either positive or negative ~\flj2 < 
X < \/7/2. The properties of the Hamiltonian of Eq. Jl6| ) 
are investigated by taking the scaled parameters 77 = 
e/[e + An{N — 1)] and x = ~^x/V7 randomly on the 
intervals < 77 < 1 and — 1 < X < 1 (these coeffi- 
cients have been used as control parameters in a study 
of phase transitions in the IBM |§,|1). In Figs. and || 
we show the corresponding probability distribution and 
correlation plot for the consistent-Q formulation of the 
IBM with realistic interactions. Although in this case 
the points are concentrated in a smaller region of the 
plot than before, the results show the same quantitative 
behavior as with random one- and two-body interactions. 
In Fig. § we have identified each of the dynamical sym- 
metries of the IBM (and the transitions between them) . 
There is a large overlap between the regions with the 
highest concentration of points in Figs. ^ and |j| 



These results, i.e. the dominance of J p = + 
ground states and the occurrence of both vibrational 
and rotational structures, are not based solely on ener- 
gies, but also involve wave function information via the 
quadrupole transitions. The use of random interactions 
(both in magnitude and sign) show that these regular 
features arise from a much wider class of Hamiltonians 
than are generally considered to be realisitic, and are, 
to a certain extent, independent of the specific character 
of the interactions. This too is in qualitative agreement 
with the empirical observation of robust features in the 
low-lying spectra of medium and heavy even-even nuclei 
This leads naturally to the question of what is 
the cause of this behavior, if the only ingredients of the 
calculations are the one- and two-body nature of the in- 
teractions, the number of bosons N and the structure of 
the model space? 

In order to see to what extent the results found above 
^lf depend on the rank of the interactions, we study 
the effect of the inclusion of three-body interactions in 
the Hamiltonian. Three-body interactions are of special 
interest in the IBM, since they can give rise to stable 
triaxial deformations |23|, which are absent in the case 
of one- and two-body interactions only. We consider the 
most general one-, two- and three-body IBM Hamiltonian 
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where Hi and H 2 are given in Eqs. (|l0|) and (|il|). H3 
contains the 17 possible three-body random interactions 
and can be written in a similar fashion. Again we find a 
dominance (61.9 %) of J p = + ground states; in 12.0 % 
of the cases the ground state has J p = 2 + , and in 17.7 
% it has the maximum value of the angular momentum 
J p = 32 + , very close to the results for the case of one- 
and two-body interactions only (63.4, 13.8 and 16.7 %, re- 
spectively). Also the probability distribution P(R) shows 
the same behavior as for one- and two-body interactions. 
In Fig. ^| we compare the result for H12 (solid curve, see 
Fig. ||) with H123. In both cases we see a very clear 
structure with pronounced peaks at the vibrational and 
rotational values of the energy ratio. Since the Hamilto- 
nian H123 depends on 26 independent random numbers 
(2 one-body, 7 two-body and 17 three-body matrix ele- 
ments), there is less correlation between its many-body 
(N = 16) matrix elements than in the case of Hi 2 . This 
results in somewhat less pronounced (lower and broader) 
- but nevertheless very clear - peaks in the probability 
distribution P{R). 



IV. SUMMARY AND CONCLUSIONS 

In this work, we discussed global properties of nuclear 
structure using random Hamiltonians. 

We first considered the problem from a shell model 
perspective, focussing on a system of identical neutrons 
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in the sd shell. We confirmed the conclusion reached 
earlier by Johnson and collaborators |ll],[l2j that nuclei 
with many nucleons favor J p = + ground states, even 
without a dominant pairing component in the force. We 
demonstrated further that this is not a consequence of 
time-reversal invariance of the random Hamiltonian. Fi- 
nally, we showed that systems of identical nucleons inter- 
acting via random two-body interactions tend to favor a 
seniority structure, in accord with conclusions reported 
recently in jl3). There is little or no evidence for the 
occurrence of vibrational and rotational bands. 

We then considered the same issues in the context of 
the IBM, a collective model that from the outset empha- 
sizes the role of monopole and quadrupole pairs. Here 
too we found that + ground states predominate, exactly 
as in the shell model analysis. In contrast, we found 
that the IBM strongly favors both vibrational and ro- 
tational structures, as evident from energy ratios of low- 
lying states and their corresponding BE(2) ratios. These 
conclusions emerged from a much wider class of Hamil- 
tonians than is usually thought to be 'realistic'. The in- 
clusion of three-body random interactions did not change 
these basic conclusions, as long as the number of bosons 
is sufficiently large. This suggests that the observed vi- 
brational and rotational features represent general and 
robust properties of the IBM model space, and do not de- 
pend significantly on details of the interaction. Since the 
structure of the model space is completely determined by 
the degrees of freedom, our results emphasize once again 
the importance of the selection of the relevant degrees of 
freedom. 

The results that we obtained with random Hamilto- 
nians in the shell model and the IBM are in qualitative 
agreement with the empirical observation of robust fea- 
tures in the low- lying spectra of medium and heavy even- 
even nuclei and their tripartite classification into senior- 
ity, anharmonic vibrator and rotor regimes @,||. The 
analysis with random interactions shows that seniority 
arises as a global property of the shell model space, while 
vibrational and rotational bands arise as general features 
of the interacting boson model space. However, the IBM 
is based on the assumption that low-lying collective ex- 
citations in nuclei can be described as a system of inter- 
acting monopole and quadrupole bosons, which in turn 
are associated with generalized pairs of like- nucleons with 
angular momentum L = and L = 2. Therefore it would 
be very important to establish whether vibrational and 
rotational features can also arise from ensembles of ran- 
dom interactions in the nuclear shell model, if appropri- 
ate (minimal) restrictions are imposed on the parameter 
space. 
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TABLE I. Percentage of ground states with angular mo- 
mentum J for GOE with N = 2 identical particles in the sd 
shell (the nucleus 18 O). 



TABLE III. Percentage of J p = + ground states for 
TBRE with N = 6 identical particles in the sd shell (the 
nucleus 22 0) as a function of e. 



N J d a Basis GOE e TBRE 



2 3 2.00 21.4 % 15.9 % o.OO 67.7 % 

1 2 1.73 14.3% 4.9% .25 69.3 % 

2 5 2.45 35.7 % 68.3 % 0.50 71.7 % 

3 2 1.73 14.3 % 6.1 % 0.75 74.0 % 

4 2 1.74 14.3 % 4.8% L00 76.8 % 



TABLE II. Percentage of ground states with angular mo- 
mentum J for TBRE with TV = 4 and 6 identical particles in 
the sd shell (the nuclei 20 ' 22 O). 



TABLE IV. Energies and B(E2) values in the dynamical 
symmetry limits of the IBM [§. In the (7(5) and SO(6) limits 
we show the result for the leading order contribution to the 
rotational spectra. 
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Basis 
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6.24 
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55.9 % 
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12 


5.05 
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21 


5.37 
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22.7 % 
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15 


4.79 


18.5 % 


1.4 % 
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15 


5.12 


18.5 % 


12.3 % 
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4.65 


7.4 % 


1.5 % 




6 
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4.69 


3.7 % 


1.3 % 



E(4+)-E(0+) B(E2;4+^2+) 
E(2+)-E(0+) B(E2;2+->0+) 



17(5) 2 2^ 

cn/fi\ 5 10(JV-l)(Af+5) 

W\IJ) 2 7JV(AT+4) 
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FIG. 1. Level distributions for N = 2 particles ( ls O). 
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FIG. 4. Probability distribution P(R) of the energy ratio 
R = [E(4+) - E(Q+)]/[E(2+) - E(Q+)] with / P(R)dR = 1 
for iV = 6 particles in the sd shell with random two-body 
interactions. 
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FIG. 2. Level distributions for N = 4 particles ( 20 O). 
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FIG. 3. Level distributions for N = 6 particles ( 22 Q). 



FIG. 5. Probability distribution P(R) of the energy ratio 
R = [£(4+) - E(Q+)]/[E(2+) - E(Q+)] with / P(R)dR = 1 in 
the IBM with random one- and two-body interactions. The 
number of bosons is N = 16. 
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FIG. 6. Correlation between ratios of B{E2) values and 
energies in the IBM with random one- and two-body interac- 
tions. The number of bosons is N = 16. 
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FIG. 7. As Fig. ^, but in the consistent-Q formulation of 
the IBM. 
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FIG. 8. As Fig. [], but in the consistent-Q formulation of 
the IBM. 




FIG. 9. As Fig. 5| but for random one- and two-body inter- 
actions _Hi2 (solid line) and random one-, two- and three-body 
interactions H123 (dashed line). 



